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Abstract 



Lattice gauge theories with Wilson fermions break chiral symmetry. In the 
U(l) axial vector current this manifests itself in the anomaly. On the other 
hand it is generally expected that the axial vector flavour mixing current 
is non-anomalous. We give a short, but strict proof of this to all orders of 
perturbation theory, and show that chiral symmetry restauration implies a 
unique multiplicative renormalization constant for the current. This constant 
is determined entirely from an irrelevant operator in the Ward identity. The 
basic ingredients going into the proof are the lattice Ward identity, charge 
conjugation symmetry and the power counting theorem. We compute the 
renormalization constant to one loop order. It is largely independent of the 
particular lattice realization of the current. 
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1 Introduction 



Any realization of fermions on the lattice has to respect the constraints imposed by 
the Nielsen- Ninomiya theorem Q. Whereas Wilson fermions break chiral symmetry 
explicitly, Ginsparg- Wilson fermions [Q, |J have an exact chiral symmetry on the 
lattice that is generated by composite local lattice operators 0. In both cases, 
the continuum chiral flavor mixing symmetry and the anomaly are to be properly 
reproduced as the cutoff is removed. 

In a recent paper |J it was shown that under very general conditions on the lattice 
Dirac operator, which, in particular, are satisfied both for Wilson and for Ginsparg- 
Wilson fermions, the axial anomaly is correctly generated in the continuum limit. 
The main conditions are gauge invariance, absense of doublers, and locality on the 
lattice in a more general sense. The origin of the anomaly is traced back to an 
irrelevant, local lattice operator in the axial vector Ward identity. 

For Ginsparg- Wilson fermions, the composite operator in the chiral transforma- 
tion which ensures an exact flavor mixing symmetry on the lattice, stays irrelevant 
under renormalization ||. As a consequence, the axial vector current does not re- 
quire renormalization. On the other hand, for Wilson fermions it is not obvious, 
although widely believed, that the chiral symmetry becomes restored in the con- 
tinuum limit. Below we give a short but strict proof of this assertion to all orders 
of perturbation theory, based on lattice power counting for massless theories 0. 
Although we explicitely refer to Wilson fermions, the result is as general as that 
for the anomaly generation given in 0. As we shall show, the only role played by 
the irrelevant, symmetry breaking operator in the flavor mixing axial vector Ward 
identity is to give rise to a unique multiplicative renormalization Zj of the axial vec- 
tor current, ensuring that chiral symmetry is restored in the continuum limit. We 
compute the one-loop contribution to Zj as a function of the Wilson parameter. The 
result is largely independent of the particular lattice regularization of the current. 

2 General framework 

Although our general proof will be given for QED, it generalizes in an obvious way 
to non-abelian gauge theories with massless fermions. 

2.1 Renormalized lattice QED 

The action for renormalized QED is given by 





x£aZ 4 fJ,^u=0 
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with g the renormalized gauge coupling constant and U{x,/i) = exp {iagA^x)) G 
U(l). The fermion action is given by 

S f = a 4 Z^(x)(D[U]+m )tlj(x), (3) 

xeaZ 4 

with if> a 2-flavor Dirac spinor field and with D[U] the Wilson Dirac operator, 

1 3 

D[U]ifj(x) = — Y^[{^-r)U(x-,tj)iP(x + aj?) 

-{in + r)U(x - dp; ^^^(x - dp) + 2rif>(x) . (4) 

m is the bare fermion mass, which for massless fields must be tuned to its critical 
value of 0(g 2 ). S g f denotes the gauge fixing action. For concreteness we choose the 
Lorentz gauge, 



S S ,(A) = a' 



E 5 (tlw) , (5) 

xGaZ 4 \/i=0 / 



with A > the gauge fixing parameter. Here and in the following, <9 M and d* denote 
the forward and backward lattice difference operators, respectively, 

%f(x) = f{x + dp) - f{x), d*J{x) = f{x) - f{x - dpi), (6) 

where // is the unit vector in \i direction. 

The generating functional W of the connected correlation functions is given by 



expW(J,r],r]) = / I difj{x)difj{x)Y\dA fJ/ {x) 

X \ [J, , 

■ exp^-S{A,if},^) + S c {A,if},^;J,r],fj)^ (7) 

with source term S c 

S c (A,iP,^;J,ri,ri)=a 4 J2{j2 J ^ A ^ + ^ ^ 

X fl 

The vertex functional T is obtained by a Legendre transformation 

W{J, 77, 77) = F{A, iP, i^) + a 4 J2 (J2 M X )M*) + v{x)Hx) + ^{x)r]{x)) , (9) 



X fl 

where 

a » / \ dW a i / \ dW a — — . . dW 

aMx) = djjx)' a ^ x) = w^r ai,{x) = ~w*)- (0) 
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By T^ n ' m \k, I) we denote the momentum space vertex function of n fermion pairs 
and m gauge fields, with their collected momenta denoted by k and I, respectively. 
For Q any composite local lattice operator, we write Fq (q; k, I) for the vertex 
function with one insertion of Q, with q its momentum. Momentum conservation is 
implied. Massless fermions require that 

trf (1 ' 0) (A; = 0) = (11) 

to be achieved by tuning mo, where the trace is taken in spinor space. Za and 
are uniquely determined by appropriate normalization conditions at non-exceptional 
momenta, e.g. by 



i d 



*dkn 



1 d 



tr l0 r^(k) r = i, -± " rS 2) (fc) r = % (12) 



k 2dk 



k 



where k = (]I ^ 0, 0, 0, 0), k = (2/a) sin(fca/2), k = (1/a) sin(fca). 
2.2 Symmetries 

Below we make explicit reference to the charge conjugation symmetry 

r(A c ,^ c ,^) = r(A,i>,i>), (is) 

where 

A°(x) = -A^x), ^ c {x) = CV(x) T , ^ C (x) = -^{xfC 1 . (14) 

The superscript T denotes transposition and C the charge conjugation matrix sat- 
isfying 

C-\C = - 7 J, /i = 0,...,3. (15) 
Furthermore, applying a gauge transformation leads to the local Ward identity 
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or , or 



-t\a d;d u d;A„(x) = 0. (16) 



It implies that the renormalized action is of the form as stated above. 

3 Chiral symmetry breaking and symmetry restora- 
tion 

Chiral symmetry is broken by the Wilson Dirac operator. Under a local, flavor 
mixing chiral transformation 

5ip(x) = ie(x) a a -f 5 ip(x), 5ip(x) = ie(x) ^(x)7 5 a a , (17) 
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where a a , a = 1,2, 3, denote the Pauli matrices acting in flavour space, the action 
transforms according to 



8S = a 



x \ n / 



18) 



with 



V a {x) = ^(x)^a a ip(x). (19) 

The gauge invariant local operators and A are not uniquely determined by (fL8|). 
In general, A is a local lattice operator which is classically irrelevant, that is, 



lim A n (x) 



0. 



a->0 



(20) 



It has UV degree 4 and IR degree 5. A convenient representation of j m and of A c 
is given by 



1 



Jfia yx) 



z i>2 (jP( x )(j» + s)^ 5 a a U(x;Li)ip(x + aft) 
+ip(x + a/x)(7 M - s)j 5 a a U(x; /i)" 1 ^) 



1 3 

A a (x) = /^2{( r ~ S ) ^(^)75CTa 

At=0 



(21) 



U (x; fJ>)ip(x + apt) + U(x — aft,; ji) ip( x ~ a V) ~ 2ip(x) 



+ (r + s) 



ip(x + aj2)U(x; jx) 1 + ip(x — ap)U(x — op,; Li) — 2ip{x) 



7 5 <7 a ^(:r) 

where r is the Wilson parameter, and s some arbitrary real but otherwise fixed con- 
stant. The following discussion on renormalization does not depend on a particular 
choice of s. 

We add to the source part of the action S c a term 



x a=l /i=0 



. X 



AJx) + 2m V a (x) 



(22) 



and denote the corresponding vertex functional by T'(A,ip,ip; G, F). Obviously, 
T'(A,ip^;G = 0,F = 0) = T(A, ip^). Then (03) implies that V satisfies the axial 
vector current Ward identity 



dV 



da A G 



fMCt 



X 



dV 
da 4 F n 



+ 



x 



r dV — dV ~\ 

I^VR^M " ^^ImW) I (23) 



+ 0(F,G). 



The functional identity (P5|) is equivalent to the infinite set of momentum space 
Ward identities 



i ^r^teM) - r<gg(M) = rt m) ( q ;k,i) + 2m r^ m \ q] k,i). (24) 

fi=0 

Here we have written TQ^(k,l) for the pure QED part, which is a linear combi- 
nation of r( n,m )(A;,Z) with 75<r a attached to the various external fermion lines, but 
with no composite operator inserted. According to the renormalization prescription 
of QED, it is UV finite and universal in the continuum limit. 

QED is already renormalized, but because of A ^ 0, the axial vector current 
requires additional renormalization. This renormalization is multiplicative. That is, 
there exists a renormalization constant Zj such that 

r&$(?;M) = z^f^kj) (25) 

is finite in the continuum limit, for all n and m. The renormalized current satisfies 
the Ward identities 

i E^5(g;A;,0 - f = M), (26) 



where 



'r& w >(g; fc, + i (z, frrj^fa k > 0] (27) 

At=0 

are the renormalized vertex functions with one A Q -insertion. Because of 

m = O(g 2 ) and Z j -l = 0(g 2 ), (28) 

the part in brackets on the right hand side of ( 57 ) is equivalently obtained by adding 
local counter terms to the lattice source action S c , that is, in (p2[), the square bracket 
is replaced by 

1 3 - 

A a (x) + 2m V a (x) + {Z j - l)-J2d; JfMa (x). (29) 

At=0 

We now show that for the particular choice of Zj, these counter terms provide 
precisely overall Taylor subtractions at zero momentum, for all correlation functions 
with one A-insertion, according to their overall ultaviolet lattice divergence degrees. 
Together with (|i2|) , because A is an irrelevant local lattice operator, this then implies 
that 

lmlt# (?;*,/) = 0, (30) 
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to all orders of perturbation theory, and for all n and m J7| . The renormalized axial 
vector current becomes conserved in the continuum limit. 

For the proof of this assertion, we recall that A is a local operator of IR degree 
5. This implies that T^'^ is continuous at zero momentum and is once con- 
tinuously differentiate at zero momentum. (These are the only vertex functions 
with one A-insertion that require overall UV subtractions, with overall (lattice) 
divergence degrees and 1, respectively.) 

First, charge conjugation symmetry implies that 

f£>) = 0. (31) 

Furthermore, for the massless theory, satisfying ([TT|) , we obtain from ( p6"D with n — 1, 
m = 

r£>) = o, (32) 

because T^'^ R (q; k, I) is at most logarithmically infrared divergent if all momenta 
are sent to zero. Again, using charge conjugation symmetry, we know that for small 
momenta 

3 

^3(95*0 = i^2l^l5°a + o(q,k) asg,fc->0, (33) 

At=0 

with (infrared) finite constant z&. Hence, order by order, Zj is uniquely determined 
by the requirement that 

r£,°i(3;A0 = <>{q,k) asg,A;-0. (34) 
This completes the proof. 



4 Axial Vector Renormalization Constant in One- 
Loop order 

To one-loop order the renormalization constant Zj is determined from the condition 



34]), where is given according to ( p7f) by 



^(1,0) 



(q; k) 



l—loop 



;(1,Q) 



(q;k) 



l—loop 



+ 2m 7 5 cr a + i(Z j -l) ^q^^^a- 



(35) 

The lattice Feynman diagrams that contribute to the first term on the right hand 
side are listed below. The vertices that correspond to the A-insertion are given in 
Appendix B. 



(36) 
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To implement condition fl3"3|) we expand the right hand side of (^) around vanishing 
momentum up to first order. After a tedious but straight forward calculation one 
finds that 

3 

/t=0 

where 7\ denotes the Taylor expansion to order 1 at zero momentum. Ca is given 
by the following (r-dependent) expression. 

with C = 1 for U(l) and C = (N 2 - 1)/(2N) for SC/(iV), and 

h a (£) = cos 4 U 2 - Il + ^ V(T (l) ) - I 2 aVa (£) 

«2 1 / » 



where 



+ 2 + 4 ^ V C ° S 2 (39) 



3 3 

I = 2sin|, £ 2 = I, = sin^, ? = ^1% 

M r = ^(/) = 2cos 2 | - ^cos 2 |. (40) 

At=0 

a is any one of the indices 0, . . . , 3. Condition (34) now determines Zj to be 

Zj = l + Ca- (41) 

Ca does not depend on a particular realization of the axial vector current j^ a , 
Eqn. (|2~ID , that is, it is independent of the parameter s. 

Note that the (finite) renormalization constant of the axial vector current is 
solely determined from diagrams involving the insertion of a classically irrelevant 
operator. The integral (|38|) is evaluated numerically. The dependence of Ca on 
the Wilson parameter is shown in Fig. |l|. From this figure we see that Ca is a 
monotonically decreasing function of r > 0. There exists no non-zero value of r for 
which the axial current is not renormalized. In particular, for the commonly used 
value r = 1 we obtain Ca = — 0.0549Cg 2 . 

5 Conclusion 

In this paper we have investigated the renormalization of the flavour mixing axial 
vector current for massless gauge theories with Wilson fermions. The corresponding 



S 



1 2 3 4 5 

r 

Figure 1: Ca(^) /(Cg 2 ) as function of the Wilson parameters r. We have Ca( — r ) — 
Ca(^) and Ca(^) 7^ whenever r ^ 0. 

axial vector Ward identity involves a symmetry breaking lattice operator A which 
is local and classically irrelvant. Using the lattice power counting theorem for mass- 
less field theories, we have shown that A uniquely determines the renormalization 
constant of the axial vector current in such a way that the chiral symmetry becomes 
restored in the continuum limit, to all orders of perturbation theory. 

We have computed the renormalization constant to one-loop order. It is largely 
independent of a particular lattice realization of the current and non- vanishing when- 
ever the Wilson parameter r ^ 0. 

Although we have considered Wilson fermions, the result of symmetry restoration 
in the continuum limit is quite general. It holds for any lattice Dirac operator that 
satisfies a general set of conditions. These conditions are gauge invariance and 
charge conjugation symmetry, absense of doublers, and locality in the more general 
sense as stated in || |3| . 

A Small momentum behavior 

The infrared properties of A ensure that the vertex functions \q;k,l) and 
\q; k) are continuous and once continuously differentiable at zero momentum, 
respectively. Their regular parts are obtained from the small momentum behavior 
of that part of the vertex functional T' that is linear in the source F. In momentum 
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space it reads 



d 4 q d 4 ki d A k 2 



+ 



W (2,). P .).5«H< 2 ')' ,( ' +t ' +fe: 



(42) 
(43) 



with c-numbers p r ,pi,r] and £. Applying the transformation QT4] ) yields the same 
expression with 

Pr^ Pi, v^v, € -» (44) 

The symmetry (|13|) thus implies that p r = and £ = 0. The vanishing of r\ is 
implied by the chiral Ward identity This implies the statements (PT|)-(P3]). 



B Feynman rules 



We state the Feynman rules for the insertion of one A Q -operator, Eqn. fl2~T|), that are 
required for the computation of the axial vector current renormalization constant 
Zj to one-loop order. For simplicity the rules are given for gauge group U(l). 





where 




[r — s) k 2 + (r + s) k\ 



J5<7 a 



(45) 
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(r - s) (q + h + k 2 )^ - (r + s) (q-h- k 2 ) fl 



(46) 



g 2 a 



(r — S) C q+ k 1 +k 2 ,fi + ( r + s ) Cq~ki-k 2 ,n 



a 

(47) 



k„a 



2 . k„a 



Ck u = cos , k u — - sin - J — 
k ,n 2 M a 2 



(48) 
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